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1 Introduction 

Structures of proton-rich nuclei are important for the rapid-proton (rp) pro- 
cess of the nucleosynthesis, which takes place in the hydrogen burning stage 
in stellar site. Since the strong interaction keeps the charge symmetry very 
well and the Coulomb energies are almost state-independent in a nuclide, 
energy spectra are quite analogous between mirror nuclei. Hence we usually 
estimate the level structures of Z > A^ exotic nuclei from their mirror part- 
ners. However, for example, the excitation energies of the 1/2^ first excited 
states in ^^C and ^^N show large discrepancy, which is called Thomas-Ehrman 
shift (TES)|1]]Q. The TES may have a significant influence on the scenario of 
the rp process, and it is highly desired to predict the TES correctly. Recent 
experiments provide us with valuable information of energy levels oi Z > N 
nuclei around ^^O. The TES has conventionally been regarded as an effect 
of the Coulomb force on a loosely bound or unbound proton occupying an 
s-orbit. With the aid of the recent data, it is being possible to argue whether 
this mechanism is sufficient to account for the TES in various mirror nuclei. 
The difference in the single-particle (s.p.) energies leads to a certain differ- 
ence in the s.p. wave functions between protons and neutrons, in general. 
This affects the matrix elements of residual nuclear interaction (RNI), even 
though the original A^A^ force is charge symmetric. The question is whether 
this effect on the TES is sizable. Since the nuclear interaction has short-range 
character, it is expected that the RNI becomes weaker as the s.p. wave func- 
tions distributes over a wider region. The RNI reduction due to the broad 
radial distribution of the s.p. wave functions typically amounts only to a few 
percent 0, which does not cause notable TES for low-lying states. However, 
because a loosely bound s-orbit can have very long tail owing to the lack 
of the centrifugal barrier, the levels involving such an s-orbit may have sub- 
stantial contribution of the RNI to the TES. In this paper we investigate the 
TES in A ~ 16 nuclei from a phenomenological viewpoint, primarily focusing 
on effects of the RNI on the TES, via the data of the i^N-^^^F, is^-isp ^^^ 
i6Q_i6]\j-g mirror pairs. 



^ There has been a confusion in the term 'Thomas-Ehrman shift'. In some references 
it was used in a restrictive meaning of a specific effect of the Coulomb force |2|. In this 
paper we use it in more general sense of the level shift between mirror nuclei. 



2 Phenomenological study of Thomas-Ehrman 
shifts around ^^O 

2.1 16^-16f 

We shall take the (Opi/2)^"^ Cg) (0^5/2151/2)"^ model space on top of the ^^O 
inert core. For neutron-rich nuclei with Z<8<N,ni = 8 — Z and 
722 = A^ — 8, and vice versa for their mirror partners. The single-particle 
(hole) energies are obtained from the data of ^^O and ^^F (^^N and ^^0)|^. 
Taking into account their mass differences from ^^0[^, we have (in MeV) 

en(04/2) = -4.144, e„(lsi/2) = -3.273, 

ep(0rf5/2) = -0.600, ep(lsi/2) = -0.105, (1) 

and 

ep(0Pv2) = 12.128, e,(0p-/^2) = 15.664. (2) 

The shift in E^(l/2+) of ^^F (i.e. Ae^"*^ = ep(lsi/2)-ep(0rf5/2) = 0.495 MeV) 
from that of ^^O (i.e. Ae^-"^ = en(lsi/2) -en(0rf5/2) = 0.871 MeV) is a typical 
TES. Because the proton in the lsi/2 orbit is loosely bound and free from the 
influence of the centrifugal barrier, its wave function spreads in a radial direc- 
tion (like halo or skin structure), leading to weaker Coulomb repulsion than 
that of 0(i5/2. This difference in the Coulomb energy gives rise to the TES for 
such a core plus one-particle system [p. The mechanism how the energy shift 
of Aef,"'^ from Ae^^'^ occurs has recently been investigated in some detail in 
Ref.p]. The observed energy spectra of ^^N and ^^F, the T^ = ±1 mirror 
nuclei, show a remarkable difference^]. Even the ground state spins do not 
match, being 2~ in ^^N and 0~ in ^^F. On top of the ^^O core, the lowest 
four states in these nuclei (2^, 0^, 3^, 1~ in ^^N and 0~, 1^, 2~, 3~ in ^^F) are 
classified into the |0pr,20(i5/2; J = 2~,3") and |0pr/2lsi/2; J = 0^, 1^) mul- 
tiplets. The difference in low-lying energy spectra is linked to the relatively 
low energy of the proton lsi/2 orbit. The smaller Aef,"'^ than Ae^~'^ in Eq. (|l]) 
shifts down the 0~ and 1~ states of ^^F. However, the amount of the observed 
TES in the 0~ and 1~ states is ~ 0.70 MeV on average, notably larger than 
Ae^-'^ - Ae^-'^ = 0.376 MeV. It is likely that the two-body RNI has a sub- 
stantial contribution to this TES. Since the last proton is unbound in ^^F 
while bound in ^^F, the relative energy of (lsi/2)p may be further lowered 



in ^^F, mainly by the effect of the Coulomb force. However, it is not sim- 
ple to evaluate separately the RNI and the nucleus- dependence of ep(lsi/2) 
(or Aep"''). We here assume the e's of Eqs. (^0), ignoring the nucleus- 
dependence of Aep"'', whose argument will be given later. Then the matrix 
elements of the residual proton-neutron interaction between a Opi/2 hole and 
an s, d particle can be derived from the experimental levels of ^^N and ^^F. 
The results are presented in Table |l]. We here denote the diagonal matrix ele- 
ment ((ji)p(j2)r; ^|^|(ji)p(j2)r; J) by Vpr{iii2] J), where (p, r) = (p,n). For 
example, Kpn (0^^/20*^5/2 ; ^) = ((0p7/2)p(0'^5/2)n; ^|V^|(0j9^/2)p(0'^5/2)n; ^) and 

Vnp{0p^^2^d5/2; J) = ((0p-/2)n(0rf5/2)p;J|^|(0p^/2)n(0rf5/2)p;J). While the 

matrix elements Vnp (0^7/30^5/2; J) deduced from ^^F are similar to V5pn( 0^7,30^5/2; J) 
from ^^N, the elements regarding lsi/2, Vpn(0p7,2lsi/2; J) and Kip(0p7/2lsi/2; J), 
show obvious discrepancy. The Vnp(0p7/2lsi/2; J) element is smaller by a fac- 
tor of about 0.7 than ypn(0p7,2lsi/2; J), both for J = 0^ and 1^. The reduc- 
tion of the proton-neutron RNI matrix elements indicates that the TES in 
the ^^N-^^F pair owes a part to the nuclear force, not only to the Coulomb 
force which relatively shifts down ep(lsi/2). The reduction of Kip compared 
with Vpn should originate in the difference of the s.p. radial wave functions 
between a proton and a neutron. The Vpp-Km difference in the RNI has been 
investigated in a wide mass region[Q], which yields a few percent reduction 
of Vpp relative to V^n, as a result of the proton-neutron difference in the 
s.p. wave functions. This coincides with the Kip/Vpn ratio for 0^5/2 shown 
in Table |I|. On the other hand, the present RNI reduction with respect to 
lsi/2 is remarkably stronger than the global systematics. The strong reduc- 
tion of the RNI is possibly an effect of the loosely bound s-orbit. Because of 
the lack of the centrifugal barrier, the radial function of the s-orbit depends 
appreciably on the separation energy. Since the proton lsi/2 orbit is loosely 
bound in the proton-rich nuclei of this mass region, the lsi/2 proton wave 
function Ru^^^^rp) distributes with a long tail and its spatial overlap with 
another nucleon is depressed, in comparison with Ris-^.^i^r^). Therefore nu- 
clear force is expected to give appreciably smaller matrix elements if loosely 
bound or unbound protons are involved. We shall examine this mechanism in 
Section ^. It is noted that the residual interaction relevant to the low-lying 
levels is repulsive in ^^N-^^F, because of the particle-hole conjugation. In 
addition to the small Ae^''^ relative to Ae^""^, the reduction of the repulsive 
RNI lowers the levels involving (lsi/2)p, relative to those having (0(i5/2)p. 



Thus the TES is enhanced in this pair of mirror nuclei. We shall investigate 
other mirror pairs based on the above phenomenological Hamiltonian. The 
validity of the current approach would be tested via the systematic study 
of this sort, whereas for the time being we put aside the possibility of the 
nucleus-dependence of e's. As will be shown, the systematic study supports 
the RNI reduction picture on the TES. 

2.2 15Q_15p 

By using the empirical e's of Eqs. (|l],|D and the (V^)'s of Table 0, we calculate 
energies of the low-lying levels 5/2+ and 1/2+ of ^^C and ^^F (T^ = ±3/2) in 
the (0pi/2)~^ ® (0(i5/2lsi/2)^ model space. The calculated energy levels are 
shown in Fig. |l] together with the experimental data[Q. The level inversion 
occurs; the 1/2+ states, instead of 5/2+, become lowest for both nuclei. This 
inversion is reproduced by the shell model Hamiltonian, due to the repulsion 
shown in Table |I| which is stronger between OpT/g and 0^5/2 than between 
0^-/2 and lsi/2. The 5/2+ level is observed at E^ = 0.740 MeV in ^^C, while 
at 1.300 MeV in ^^F. The shell model yields ^^(5/2+) = 0.563 MeV for ^^C 
and 1.400 MeV for ^^F. The TES in the i^C-^^F pair is thus described with 
a reasonable accuracy, though slightly overshot, within the framework of the 
phenomenological shell model. Weaker repulsion in \4p(0p^/2lsi/2; J) than 
in VpniOPih^si/2; J) plays an appreciable role in the TES. The Vpp{0p^?2J J = 
0+) and Kin(0pr/2i ^ ~ ^^) Kiatrix elements, which do not affect the excita- 
tion energies, can be evaluated from the ground-state energies of ^^C and ^^O. 
We can then calculate the absolute values of the energies, not only the exci- 
tation energies, of the ^^C and ^^F levels. The biggest discrepancy is found 
in ii^(l/2+) of ^^C, which is overestimated by 0.166 MeV, whereas the other 
energies are reproduced within the 0.1 MeV accuracy. This may suggest that 
an additional effect is missed for the lsi/2 neutron, whose separation energy 
is small (1.218 MeV) in ^^C. 

2.3 i^C-i^Ne 

The ^^C-^^Ne pair (T^ = ±2) is significant as well, in investigating the ef- 
fect of the RNI on the TES. The low-lying states of these nuclei have the 
0p7/^2 ® (0(i5/2lsi/2)^ configuration. Since the 0^7/^2 P^^^ does not contribute 



to the excitation energy, the TES can disclose difference between proton- 
proton (VJpp) and neutron-neutron (V^n) interactions in the sii-shell. As an 
effective interaction in the sd-shell, the so-called USD interaction has suc- 
cessfully been used. Although the USD interaction is derived for the full 
sd-shell calculation, we neglect the 0^3/2 components, since the 0^3/2 orbit 
is hardly occupied in low- lying states of the nuclei around ^^O. Indeed, we 
can well reproduce the low-lying levels of ^^O with the USD interaction in 
the {0d5/2lsi/2Y space. We carry out the shell model calculation in the 
0p^?2 ® (0<i5/2l'Si/2)^ space, with the Hamiltonian comprised of the empirical 
e's and (V^)'s (see Eqs. (|l|,0) and Table |l|) as well as of the USD matrix ele- 
ments. The binding energy of ^^C is reproduced with the accuracy of about 
0.1 MeV. In computing the binding energy of ^^Ne, the residual two-body 
Coulomb interaction has to be estimated. With the s.p. wave functions 
in the Woods-Saxon potential which will be mentioned in Section ^, the 
two-body Coulomb force yields approximately constant energy shift of about 
0.4 MeV for low-lying levels, within the accuracy of 0.1 MeV. If we use the 
charge-symmetric (i.e. l^p = l^nn) USD interaction with this Coulomb cor- 
rection, the binding energy of ^^Ne is seriously overestimated by as much 
as 0.8 MeV. Because of the level inversion in ^^C and ^^F, lsi/2 lies lower 
than 0(i5/2 in an effective sense. Thereby the ground state consists mainly of 
the Op7/2 ® l^irt configuration, with small admixture of OpT/g (S> Orfg^. It is 
likely for the RNI matrix elements involving (lsi/2)p to suffer some amount 
of reduction, because the lsi/2 protons are bound loosely (or unbound). For 
this reason we reduce the USD matrix elements concerning the (lsi/2)p orbit 
by an overall factor C,s, while not changing the other matrix elements. The 
binding energy of ^^Ne is found to be reproduced if we set C,s — 0.7 (i.e. 
Vpp{0d5/2lsi/2; J) — 0.7 X Kin(0(i5/2lsi/2; J), and so forth). It is notable that 
this factor is close to the proton-neutron ratio Kip/Vpn concerning lsi/2 in 
Table |l|. Recently ^'^^(O^) of ^^Ne has been reported^], indicating a large 
TES. ^^(0^) of ^^Ne is lower by about 1 MeV than that of ^^C. In Fig. 
the results of the ^s = 1 (i-e. no reduction of the RNI) and ^^ = 0.7 cases 
are compared with the experimental data. As has been noticed, the lsi/2 en- 
ergy relative to £(0^5/2) is deeper for protons than for neutrons. This tends 
to lower the ground state of ^^Ne, whose main configuration is 0p^?2 ® l'Si/2- 
Thus, if we use the charge-symmetric USD interaction (i.e. ^^ = 1), E^{02) of 
^^Ne becomes necessarily higher than that of ^^C. The recent data of £'a;(0^) 



clearly favors the reduction of the RNI regarding the (lsi/2)p orbit. The 
reduction with C,s = 0.7 almost reproduces £'^(02') of ^^Ne, as well as the 
binding energy. Note that the residual Coulomb force does not contribute 
seriously to the TES, as far as the number of the valence protons is not large. 
With this reduction of C,s = 0.7 the known energy spectra of the ^^Ne0 are 
also reproduced. 

2.4 Discussion — Ae*"^ decrease vs. RNI reduction 

We now consider the possibility of the nucleus-dependence of Ae*""^ = ep(lsi/2)- 
ep(0(i5/2)- In extracting the RNI matrix elements from ^^N-^^F, we have as- 
sumed the s.p. energies taken from the ^^0-^^F data. One may argue that 
in the ^^F nucleus (lsi/2)p receives less Coulomb repulsion than in ^^F, be- 
cause the last proton is unbound, and that the TES in ^^N-^^F should be 
ascribed to the corresponding lowering of ep(lsi/2) relative to Cp (0^5/2), in- 
stead of the reduction of the RNI. As far as we view only the ^^N-^^F data, 
the RNI reduction does not seem to be an exclusive solution. In this regard, 
the TES in the ^^C-^^Ne pair has particular importance. Since ^^Ne has 
negative S2p (two proton separation energy), ep(lsi/2) (relative to ep(0c?5/2)) 
goes down in ^^Ne from the value obtained in ^^F, if we follow the argument 
of the nucleus-dependence of Aef,"'^ as in ^^F. However, this makes £^3,(02") 
in ^^Ne even higher than in the ^s = 1 case of Fig. |^, where we already have 
the wrong sign of the TES. It is noted that the RNI derived from ^^N-^^F 
(shown in Table |I|) is repulsive because of the particle-hole nature, while the 
sd-shell interaction crucial to the TES in the ^^C-^^Ne pair is attractive. 
Thereby, the Aep"'^ decrease and the RNI reduction gives opposite effects 
on the loosely bound (or unbound) s-orbit in ^^C-^^Ne. If only the Ae!,""^ 
decrease is considered, obvious contradiction to the data results. Thus, lower 
Erc{02) in ^^Ne than in ^^C cannot be described without the reduction of the 
RNI for (lsi/2)p. The data implies that, although the nucleus- dependence 
of Aep"'' may exist, its effect on the TES seems much less significant than 
that of the RNI reduction. On the contrary, the reduction of the RNI natu- 
rally accounts for the TES's around ^^O, in particular those of ^^N-^^F and 
i6Q_i6j^g^ simultaneously. 



2.5 ^^,(3+) of i^Ne 

The TES is not apparent for the lowest-lying states of the ^^0-^^Ne mirror 
pair, since their main configuration is Orf^/g. On the other hand, the lowest 
3^ state, which is observed at E^ = 5.378 MeV in ^^O, is expected to have 
the 0(i5/2lsi/2 configuration. Because the lsi/2 orbit is relevant, the TES 
for this state may be sizable. The energy of the 3"'" state of ^^Ne is impor- 
tant in the scenario of the rp process[10|. Whereas an earlier experiment 



gives Ex = 4.56 MeV[|rT|, this 3+ state has not been established by recent 
experiments jT^ . By using the USD interaction with the reduction factor 
$,s = 0.7 for (lsi/2)p (together with ep's deduced from ^^F), we evaluate TES 
of 0.86 MeV for this 3+ state. This leads to ^^(3^) ~ 4.5 MeV, in good 



agreement with Ref . [pT 



3 Mechanism of RNI reduction 

We next study the mechanism of the RNI reduction concerning the proton 
lsi/2 orbit, from a qualitative (or semi-quantitative) standpoint. As has been 
pointed out, it is likely that the RNI reduction is connected to the broad dis- 
tribution of Ris^.^i^Tp). The amount of the reduction, however, is notably 
large, compared with the global trend which has been estimated to be a 
few percent 0. It is a question whether Ris-^^^{rp) distributes so broadly as 
to give RNI reduction by about 30%, despite the presence of the Coulomb 
barrier. It is not an easy task to estimate microscopically the RNI matrix ele- 
ments to a good precision. Instead, we view proton-neutron ratio of the RNI 
matrix elements (Vnp/^n and V^p/Kin) of the M3Y interaction. The M3Y 
force[|l^ basically represents the G-matrix and therefore somewhat realistic, 
and enables us to avoid tedious computation. To take into account effects 
of the loose binding with the centrifugal barrier and the Coulomb barrier, 
the single-particle wave functions are obtained under the Woods-Saxon (WS) 
plus Coulomb potential. We adopt the WS parameters of Ref. [jl4[ at ^^O, 
varying the WS potential depth Vq around the normal value —51 MeV. Even 
if the absolute values of the RNI matrix elements are not reliable, the proton- 
neutron ratios carries certain information, because they depend mainly on 
the proton-neutron difference of the s.p. wave functions. Note that core 
polarization effects are not taken into consideration in this WS-I-M3Y calcu- 



lation, which should be included in the shell model interaction. The present 
proton-neutron ratios give only qualitative (or semi-quantitative) nature of 
the RNI, since they do not need to match precisely the shell model ones (for 
instance, Table |l]). The WS potential with Vq = —53 gives (in MeV) 

en(04/2) = -7.52, e„(lsi/2) = -4.80, 

ep(0rf5/2) = -3.90, ep(lsi/2) = -1.46. (3) 

As the potential becomes shallower, the s.p. orbits are bound more loosely. 
Indeed, at Vq = —51 we have 

611(04/2) = -6.36, en(lsi/2) = -3.98, 

ep(0rf5/2) = -2.80, ep(lsi/2) = -0.76, (4) 

and at Vq = -49 

e„(04/2) = -5.23, e„(lsi/2) = -3.22, 

ep(04/2) = -1.75, ep(lsi/2) = -0.14. (5) 

Whereas the wave function is insensitive to e for the deeply bound orbits, 
it is not the case for the loosely bound orbit (lsi/2)p- The variation of 
the wave function is typically measured by the mean radius of the orbit 
TpiJ) = \l{{j)p\r'^\{j)p) (p = P,n). By varying the WS parameter Vq, we 
see how the RNI as well as rp{j) behave as e changes. For the M3Y matrix 
elements, the proton-neutron ratios Kp/Vpn with respect to the (0pi/2)~^ ® 
(0(i5/2lsi/2)^ two-body states and \4>p/Vnn with respect to the (0(i5/2lsi/2)^ 
states are depicted in Fig. 0. Though the ratios of the off-diagonal ele- 
ments {{lsl,^)p-Q+\Vmi,^)p-Q+) and m^,2lsy2)p,'2+\Vmi^)p-2+) are 
not shown, they behave in a similar manner to those of diagonal elements 
involving lsi/2. The proton-neutron ratios of the rms radii of the s.p. or- 
bits are also presented in Fig. ^ for j = 0^5/2 and lsi/2. As is expected, 
i?isj,2(rp) distributes over a broader region than Ris_^ (r^), to a certain ex- 
tent. In Fig. 1^, the rms radius of (lsi/2)p is larger by about 10-20% than 
that of (lsi/2)n, somewhat depending on Vq. In contrast, the rms radius of 
(0(i5/2)p differ only by a few percent from that of (0^5/2)11, insensitive to Vq. 
From Fig. ^ we confirm the following two points: (a) the RNI reduction well 
correlates to the increase of the rms radii of the relevant orbits, and (b) the 



matrix elements involving (lsi/2)p can be reduced from those of (lsi/2)n by as 
much as a few tens percent around ^^O. The former point is consistent with 
Ref ■ |]T5| , though we use more realistic s.p. wave functions (but less realistic 



G-matrix) than in Ref.|jT5l|. The latter implies that the broad distribution 
of i?isj^2(^p) seems accountable for the RNI reduction. Although there re- 
main additional interests in the RNI reduction (e.g. accurate estimate of the 
reduction factor, nucleus- and/or state-dependence of the reduction factor), 
they will require reliable treatment of the core polarization effects, which is 
beyond the scope of the current study. We just point out at this moment 
that, due to the broad distribution of the s.p. wave function, the core polar- 



ization effects tend to diminish ||15|| and therefore the shell model interaction 
may be reduced further. The residual Coulomb force hardly contributes to 
the excitation energies of low-lying states, for the nuclei around ^^O, as far as 
the number of valence protons remains small. The state-dependence of the 
residual Coulomb force is less than 0.1 MeV for the nuclei under discussion, 
if estimated with the above WS wave functions. 



4 Summary 

The Thomas- Ehrman shifts generally occur in the A ~ 16 region, where 
the lsi/2 proton is unbound or loosely bound. As well as the difference 
between Ae^~'^ and Aef,"*^, the reduction of the residual nuclear interaction 
matrix elements involving the lsi/2 proton plays an important role in the 
TES. As has been deduced from the nuclei ^^N and ^^F, the matrix elements 
V^p(0pr,2lsi/2) is notably smaller than Vpn{0p^X^si/2), by a factor of about 
0.7. This factor is remarkably smaller than the general trend of the proton- 
neutron asymmetry in the RNI. Similar reduction of Vpp in the sd-shell (rel- 
ative to Van) accounts for the TES in £'a;(0^) of the ^^C~^^Ne pair as well 
as the mass of ^^Ne. It is remarked that the RNI reduction is far more sig- 
nificant than the nucleus-dependence of Aep~'^, as is argued in connection to 
ExiOt) of ^^Ne. Taking into account the RNI reduction, the TES's observed 
in i5Q_i5p g^j^^ other pairs are understood within the phenomenological shell 
model. On the same ground the astrophysically important E^iSf) of ^^Ne is 
predicted to be ~ 4.5 MeV. The reduction of the residual interaction seems 
to originate in the broad radial distribution of wave function of the lsi/2 
proton, which is bound loosely (or unbound) and is not affected by the cen- 



trifugal barrier. This picture is supported by viewing the proton-neutron 
ratio of the M3Y interaction matrix elements, which are evaluated with the 
single-particle wave functions under the Woods-Saxon plus Coulomb poten- 
tial. 

Discussions with S. Kubono, K. Kato and S. Aoyama are gratefully acknowl- 
edged. 
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Table 1: Matrix elements of residual proton-neutron interaction Vpa{jiJ2', J) 
and Kip(jij2; J) deduced from ^^N and ^^F (MeV), and their ratio. 
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Figure 1: Experimental and calculated energy spectra of the i5Q_i5p j^iirror 
nuclei. 
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Figure 2: Experimental and calculated (with and without the reduction fac- 
tor ^s) energy spectra of the ^^C-^^Ne mirror nuclei. 
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Figure 3: Proton-neutron ratios of the RNI diagonal matrix elements in the 
WS+M3Y model, for the WS potential depth varied around Vq = —51 MeV. 
The ratios involving the lsi/2 orbit are shown at the right panel, and those 
without lsi/2 but with 0^5/2 are at the left panel. The J values of the 
two-body states are indicated in the graph. The corresponding ratios with 
different Vq are connected by thin lines. The proton-neutron ratios of the 
rms radii of the s.p. orbits are also shown for j = Od^/2 (left panel) and lsi/2 
(right panel), by filled diamonds linked by thick dashed lines. 
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